Comparative functional genomics and the bovine macrophage response to strains of the Mycobacterium genus

Mycobacterial infections are major causes of morbidity and mortality in cattle and are also potential zoonotic agents with implications for human health. Despite the implementation of comprehensive animal surveillance programs, many mycobacterial diseases have remained recalcitrant to eradication in several industrialized countries. Two major mycobacterial pathogens of cattle are Mycobacterium bovis and Mycobacterium avium subspecies paratuberculosis (MAP), the causative agents of bovine tuberculosis (BTB) and Johne's disease (JD), respectively. BTB is a chronic, granulomatous disease of the respiratory tract that is spread via aerosol transmission, while JD is a chronic granulomatous disease of the intestines that is transmitted via the fecal-oral route. Although these diseases exhibit differential tissue tropism and distinct complex etiologies, both M. bovis and MAP infect, reside, and replicate in host macrophages - the key host innate immune cell that encounters mycobacterial pathogens after initial exposure and mediates the subsequent immune response. The persistence of M. bovis and MAP in macrophages relies on a diverse series of immunomodulatory mechanisms, including the inhibition of phagosome maturation and apoptosis, generation of cytokine-induced necrosis enabling dissemination of infection through the host, local pathology, and ultimately shedding of the pathogen. Here, we review the bovine macrophage response to infection with M. bovis and MAP. In particular, we describe how recent advances in functional genomics are shedding light on the host macrophage-pathogen interactions that underlie different mycobacterial diseases. To illustrate this, we present new analyses of previously published bovine macrophage transcriptomics data following in vitro infection with virulent M. bovis, the attenuated vaccine strain M. bovis BCG, and MAP, and discuss our findings with respect to the differing etiologies of BTB and JD

GOexpress: an R/Bioconductor package for the identification and visualisation of robust gene ontology signatures through supervised learning of gene expression data

Background: Identification of gene expression profiles that differentiate experimental groups is critical for discovery and analysis of key molecular pathways and also for selection of robust diagnostic or prognostic biomarkers. While integration of differential expression statistics has been used to refine gene set enrichment analyses, such approaches are typically limited to single gene lists resulting from simple two-group comparisons or time-series analyses. In contrast, functional class scoring and machine learning approaches provide powerful alternative methods to leverage molecular measurements for pathway analyses, and to compare continuous and multi-level categorical factors. Results: We introduce GOexpress, a software package for scoring and summarising the capacity of gene ontology features to simultaneously classify samples from multiple experimental groups. GOexpress integrates normalised gene expression data (e.g., from microarray and RNA-seq experiments) and phenotypic information of individual samples with gene ontology annotations to derive a ranking of genes and gene ontology terms using a supervised learning approach. The default random forest algorithm allows interactions between all experimental factors, and competitive scoring of expressed genes to evaluate their relative importance in classifying predefined groups of samples. Conclusions: GOexpress enables rapid identification and visualisation of ontology-related gene panels that robustly classify groups of samples and supports both categorical (e.g., infection status, treatment) and continuous (e.g., time-series, drug concentrations) experimental factors. The use of standard Bioconductor extension packages and publicly available gene ontology annotations facilitates straightforward integration of GOexpress within existing computational biology pipelines.Department of Agriculture, Food and the MarineEuropean Commission - Seventh Framework Programme (FP7)Science Foundation IrelandUniversity College Dubli

RNA-seq transcriptional profiling of peripheral blood leukocytes from cattle infected with Mycobacterium bovis

Bovine tuberculosis, caused by infection with Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including gene expression microarrays and RNA-sequencing (RNA-seq), has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analyzed the peripheral blood leukocyte (PBL) transcriptome of eight natural M. bovis-infected and eight age- and sex-matched non-infected control Holstein-Friesian animals using RNA-seq. In addition, we compared gene expression profiles generated using RNA-seq with those previously generated using the high-density Affymetrix(®) GeneChip(®) Bovine Genome Array platform from the same PBL-extracted RNA. A total of 3,250 differentially expressed (DE) annotated genes were detected in the M. bovis-infected samples relative to the controls (adjusted P-value ≤0.05), with the number of genes displaying decreased relative expression (1,671) exceeding those with increased relative expression (1,579). Ingenuity(®) Systems Pathway Analysis (IPA) of all DE genes revealed enrichment for genes with immune function. Notably, transcriptional suppression was observed among several of the top-ranking canonical pathways including Leukocyte Extravasation Signaling. Comparative platform analysis demonstrated that RNA-seq detected a larger number of annotated DE genes (3,250) relative to the microarray (1,398), of which 917 genes were common to both technologies and displayed the same direction of expression. Finally, we show that RNA-seq had an increased dynamic range compared to the microarray for estimating differential gene expression

Whole-transcriptome, high-throughput RNA sequence analysis of the bovine macrophage response to Mycobacterium bovis infection in vitro

BACKGROUND: Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same MDM-extracted RNA. RESULTS: A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). Of the differentially expressed genes, 694 were common to both the sense and antisense data sets, with the direction of expression (i.e. up- or downregulation) positively correlated for 693 genes and negatively correlated for the remaining gene. Gene ontology analysis of the differentially expressed genes revealed an enrichment of immune, apoptotic and cell signalling genes. Notably, the number of differentially expressed genes identified from RNA-seq sense strand analysis was greater than the number of differentially expressed genes detected from microarray analysis (2,584 genes versus 2,015 genes). Furthermore, our data reveal a greater dynamic range in the detection and quantification of gene transcripts for RNA-seq compared to microarray technology. CONCLUSIONS: This study highlights the value of RNA-seq in identifying novel immunomodulatory mechanisms that underlie host-mycobacterial pathogen interactions during infection, including possible complex post-transcriptional regulation of host gene expression involving antisense RNA

Characterization of the Prophage Repertoire of African Salmonella Typhimurium ST313 Reveals High Levels of Spontaneous Induction of Novel Phage BTP1

In the past 30 years,Salmonella bloodstream infections have become a significant health problem in sub-Saharan Africa and are responsible for the deaths of anestimated 390,000 people each year. The disease is predominantly caused by a recently described sequence type of SalmonellaTyphimurium: ST313, which has a distinctive set of prophage sequences. We have thoroughly characterized the ST313-associated prophages both genetically and experimentally. ST313 representative strain D23580 contains five full-length prophages: BTP1, Gifsy-2D23580, ST64BD23580, Gifsy-1D23580,and BTP5. We show that commonS.Typhimurium prophages Gifsy-2, Gifsy-1, andST64B are inactivated in ST313 by mutations. Prophage BTP1 was found to be a functional novel phage, and the first isolate of the proposed new species “Salmonellavirus BTP1”, belonging to the P22virusgenus. Surprisingly,∼109BTP1 virus particlesperml were detected in the supernatant of non-induced, stationary-phase culturesof strain D23580, representing the highest spontaneously induced phage titer so farreported for a bacterial prophage. High spontaneous induction is shown to be anintrinsic property of prophage BTP1, and indicates the phage-mediated lysis of around0.2% of the lysogenic population. The fact that BTP1 is highly conserved in ST313 poses interesting questions about the potential fitness costs and benefits of novel prophagesin epidemicS.Typhimurium ST313

Fingertip force control during bimanual object lifting in hemiplegic cerebral palsy

In the present study we examined unimanual and bimanual fingertip force control during grasping in children with hemiplegic cerebral palsy (CP). Participants lifted, transported and released an object with one hand or both hands together in order to examine the effect on fingertip force control for each hand separately and to determine whether any benefit exists for the affected hand when it performed the task concurrently with the less-affected hand. Seven children with hemiplegic CP performed the task while their movement and fingertip force control were measured. In the bimanual conditions, the weight of the instrumented objects was equal or unequal. The durations of the all temporal phases for the less-affected hand were prolonged during bimanual control compared to unimanual control. We observed close synchrony of both hands when the task was performed with both hands, despite large differences in duration between both hands when they performed separately. There was a marginal benefit for two of the five force related variables for the affected hand (grip force at onset of load force, and peak grip force) when it transported the object simultaneously with the less-affected hand. Collectively, these results corroborate earlier findings of reaching studies that showed slowing down of the less-affected hand when it moved together with the affected hand. A new finding that extends these studies is that bimanual tasks may have the potential to facilitate force control of the affected hand. The implications of these findings for recent rehabilitative therapies in children with CP that make use of bimanual training are discussed

Hsp90 governs dispersion and drug resistance of fungal biofilms

Fungal biofilms are a major cause of human mortality and are recalcitrant to most treatments due to intrinsic drug resistance. These complex communities of multiple cell types form on indwelling medical devices and their eradication often requires surgical removal of infected devices. Here we implicate the molecular chaperone Hsp90 as a key regulator of biofilm dispersion and drug resistance. We previously established that in the leading human fungal pathogen, Candida albicans, Hsp90 enables the emergence and maintenance of drug resistance in planktonic conditions by stabilizing the protein phosphatase calcineurin and MAPK Mkc1. Hsp90 also regulates temperature-dependent C. albicans morphogenesis through repression of cAMP-PKA signalling. Here we demonstrate that genetic depletion of Hsp90 reduced C. albicans biofilm growth and maturation in vitro and impaired dispersal of biofilm cells. Further, compromising Hsp90 function in vitro abrogated resistance of C. albicans biofilms to the most widely deployed class of antifungal drugs, the azoles. Depletion of Hsp90 led to reduction of calcineurin and Mkc1 in planktonic but not biofilm conditions, suggesting that Hsp90 regulates drug resistance through different mechanisms in these distinct cellular states. Reduction of Hsp90 levels led to a marked decrease in matrix glucan levels, providing a compelling mechanism through which Hsp90 might regulate biofilm azole resistance. Impairment of Hsp90 function genetically or pharmacologically transformed fluconazole from ineffectual to highly effective in eradicating biofilms in a rat venous catheter infection model. Finally, inhibition of Hsp90 reduced resistance of biofilms of the most lethal mould, Aspergillus fumigatus, to the newest class of antifungals to reach the clinic, the echinocandins. Thus, we establish a novel mechanism regulating biofilm drug resistance and dispersion and that targeting Hsp90 provides a much-needed strategy for improving clinical outcome in the treatment of biofilm infections

Effects of external nutrient sources and extreme weather events on the nutrient budget of a Southern European coastal lagoon

The seasonal and annual nitrogen (N), phosphorus (P), and carbon (C) budgets of the mesotidal Ria Formosa lagoon, southern Portugal, were estimated to reveal the main inputs and outputs, the seasonal patterns, and how they may influence the ecological functioning of the system. The effects of extreme weather events such as long-lasting strong winds causing upwelling and strong rainfall were assessed. External nutrient inputs were quantified; ocean exchange was assessed in 24-h sampling campaigns, and final calculations were made using a hydrodynamic model of the lagoon. Rain and stream inputs were the main freshwater sources to the lagoon. However, wastewater treatment plant and groundwater discharges dominated nutrient input, together accounting for 98, 96, and 88 % of total C, N, and P input, respectively. Organic matter and nutrients were continuously exported to the ocean. This pattern was reversed following extreme events, such as strong winds in early summer that caused upwelling and after a period of heavy rainfall in late autumn. A principal component analysis (PCA) revealed that ammonium and organic N and C exchange were positively associated with temperature as opposed to pH and nitrate. These variables reflected mostly the benthic lagoon metabolism, whereas particulate P exchange was correlated to Chl a, indicating that this was more related to phytoplankton dynamics. The increase of stochastic events, as expected in climate change scenarios, may have strong effects on the ecological functioning of coastal lagoons, altering the C and nutrient budgets.Portuguese Science and Technology Foundation (FCT) [POCI/MAR/58427/2004, PPCDT/MAR/58427/2004]; Portuguese Science and Technology Foundation (FCT

Computing Power and Sample Size for Case-Control Association Studies with Copy Number Polymorphism: Application of Mixture-Based Likelihood Ratio Test

Recent studies suggest that copy number polymorphisms (CNPs) may play an important role in disease susceptibility and onset. Currently, the detection of CNPs mainly depends on microarray technology. For case-control studies, conventionally, subjects are assigned to a specific CNP category based on the continuous quantitative measure produced by microarray experiments, and cases and controls are then compared using a chi-square test of independence. The purpose of this work is to specify the likelihood ratio test statistic (LRTS) for case-control sampling design based on the underlying continuous quantitative measurement, and to assess its power and relative efficiency (as compared to the chi-square test of independence on CNP counts). The sample size and power formulas of both methods are given. For the latter, the CNPs are classified using the Bayesian classification rule. The LRTS is more powerful than this chi-square test for the alternatives considered, especially alternatives in which the at-risk CNP categories have low frequencies. An example of the application of the LRTS is given for a comparison of CNP distributions in individuals of Caucasian or Taiwanese ethnicity, where the LRTS appears to be more powerful than the chi-square test, possibly due to misclassification of the most common CNP category into a less common category